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Abstract 



In this report we present simulated event numbers, for various MC generators, for 
pion production in CC reactions on 16 0. For the simulation we used four different 
neutrino interaction generators: GENIE, FLUKA, NEUT, and NuWro, as proposed 
during the 45th Karpacz Winter School on neutrino interactions [1]. First we give a 
brief outline of the theoretical models relevant to pion production. We then present 
results, in the form of tables showing the occupancy of primary and final state pion 
topologies, for all the generated samples. Finally we compare the results from the 
different generators and draw conclusions about the similarities and differences. For 
some of the generators we explore the effect of varying the axial mass parameter or 
the use of a different nuclear model. 
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1 Introduction 

Understanding Charged-Current (CC) neutrino-nucleus interactions in the few GeV re- 
gion is very important for many current and future neutrino experiments. The study of 
neutrino-nucleus reactions in this region is complicated and requires many intermediate 
steps, such as a description of the nuclear model, understanding the neutrino-nucleon 
cross sections, modeling of hadronization, as well as the modeling of intranuclear hadron 
transport and other secondary interactions. These can all play a significant role in how we 
understand the nature of neutrinos as well as providing useful information about nuclear 
phenomena. Because of this there are a number of Monte Carlo generators and numerical 
packages dedicated to the description of neutrino interactions: GENIE [2], GiBUU [3], 
FLUKA [l], NEUT [5], NuWro [ ] and Nuance [ ] represent a large fraction of such 
generators and were all presented at the Ladek Winter School [1]. 

During the school we undertook a project to compare the predictions made by different 
generators. Of the generators mentioned previously we looked at GENIE, NEUT, Nu- 
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ance 1 , FLUKA and NuWro. For now GiBUU is not included in this study because it was 
difficult to find a consistent way of comparing it to the other simulation packages 2 . It was 
decided that we would focus our investigations on the production of pions in neutrino- 
nucleus interactions 3 . This was because they form an important background in many 
neutrino oscillation experiments and are also theoretically challenging, due to processes 
such as final state interactions (FSIs). A special web site of "Ladek MC Project" was set 
up to collate and discuss the results [8]. 

For each generator we produced a sample of mono-energetic (1 GeV) interactions 
on 16 O, using the default settings for each generator. To simplify the analysis we con- 
sider only charged current (CC) interactions. Each generator had the following processes 
enabled: quasi-elastic (QE) scattering, resonance (RES) production, deep-inelastic scat- 
tering (DIS) and coherent (COH) pion production. We then analyzed the samples for 
each generator by looking at the various pion topologies before and after any secondary 
interactions. 

Before presenting the results of the study we give a brief outline of the theoretical models 
of relevance to CC neutrino induced pion production. 



2 Overview of Neutrino-induced Pion Production 

As mentioned previously, the modeling of neutrino-nucleus interactions is complex and 
requires linking together many different pieces of theory. Here we focus on neutrino- 
nucleon cross-sections and how they are embedded in a nuclear environment. We leave 
the description of the hadronization and final state interaction models, often specific to 
a particular generator, to chapter 3. 

The total cross section for neutrino-nucleon scattering has the following form [ 

a l°N = a vN © a uN © ° 2 vN © • • • © ^In © • • • © a uN S ■ W 

The region of neutrino energies around 1 GeV is particularly troublesome. It is in this 
region that many of the above cross sections are similar in magnitude. Here resonance 
single pion production contributes ~ 30% to the total cross section [10], similar to the 
contributions from QEL and DIS processes. This is a problem experimentally as RES 
events can have indistinguishable signatures to DIS events in a detector, making it hard 
to measure each process exclusively. 

1 The results of the simulation using Nuance alongside a description of the generator are pre- 
sented in the appendix, Table 10. 

2 GiBUU uses a more sophisticated model in which the propagation/collision of resonances is 
handled explicitly. In this model the primary hadronic system is the QE nucleon, the resonances, 
and any pions from the non-resonant background. And so comparison with the other simulation 
packages is difficult as they decay the resonances before rescattering and so define the initial 
state as the decay products of all resonances. 

3 The production of other mesons, like rj or p, is also possible but in general pion production 
dominates. 



3 



Charged current QEL scattering of a neutrino on a free nucleon (u£+N— >£+N') is usually 
described using the Llewellyn Smith formalism [ ]. Although no pions are produced 
directly it is possible, through FSIs, to produce them in the final state system. Inside the 
nucleus, hadrons can be scattered elastically or inelastically, can be absorbed or charge 
exchanged and even produce extra pions (pion production). Thus, a small number of 
events with pions in the final state are expected from CCQE events, even though no 
pions were produced initially. 

The dominant CC processes that produce pions directly are DIS, COH and RES produc- 
tion. These are shown in Fig. 1. 




^ [RES] 




Fig. 1. The main types of charged current muon neutrino scattering on a free nucleon/nucleus 
that produce pions directly. From top left to bottom right are: Deep Inelastic Scattering (DIS), 
Coherent pion production, and Resonance production (RES). In the figure N is a nucleon, A is 
a nucleus and X represents the hadronic system excluding pions. 

As first proposed by Bodek and Ritchie [12] structure functions are used to describe DIS. 
Recently some progress in this field has been made using the "higher twist" QCD tech- 
nique [13]. Neutrinos can also interact with the whole nucleus (instead of with individual 
nucleons as in the previous two processes) coherently producing pions. Typically COH 
pion production is described using the original Rein and Sehgal model [11] with updates 
taking into account lepton mass terms [J ]. A further description of COH scattering is 
presented in [16] and [17]. Resonant events are usually described using the Rein-Sehgal 
model [ 8] describing the excitation of baryon resonances and pion production. 

So far we have listed cross-section models which describe the scattering of neutrinos 
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off free nucleons 4 . It is necessary to take into account the fact that these nucleons are 
not free but rather exist as bound states within a nuclear environment. The common 
approach within MC generators is to use the Relativistic Fermi gas (RFG) model where 
the Fermi motion of individual nucleons is taken into account. However its implementation 
often differs for different neutrino generators. In several papers [19] the importance of 
considering Pauli blocking and FSI effects for v e + 16 O — > eX reactions (also applying 
to f M 16 0) are shown. Also, it is shown by 0. Benhar et al. that the RFG model does 
not agree well with experimental data. A better description is offered through the use of 
spectral functions [20], as measured in electron scattering experiments. 

It is also necessary to describe hadronization, as well as the propagation of secondary 
particles out of the nucleus. The simulation must cover a description of both rescattering 
and absorption effects. A report on the modeling of final state interactions and the use of 
intranuclear cascade models was presented at the school [21]. These are often individual 
features of a generator and are described in the next section. 



3 Results of Simulations for the Different Generators 

For each of the generators we produced similar sets of 500,000 events simulating 1 GeV 
Up charged-current interactions on 16 0. We now present the results of these simulations 
in the form of tables showing the occupancy of primary and final state pion topologies. 
A primary state is defined as the topology of particles produced by the primary neu- 
trino interaction and the final state is defined as the topology of the particles after any 
secondary interactions, such as intranuclear rescattering, have taken place. We do this 
separately for all of the generators. To aid fair comparison between the tables they are 
each preceded by a brief outline of the models and physics choices in each generator 5 . 



3.1 GENIE 

GENIE simulates neutrino interactions, for all neutrino flavors and all nuclear targets, 
over a large energy range from a few MeV to several hundred GeV. The physics models 
used can, broadly speaking, be split into models which describe cross-sections, hadroniza- 
tion, and nuclear physics. Full information on all the models and physics choices used in 
GENIE can be found at [2]. 

For cross-sections: Charged current quasi-elastic scattering is modeled using an implemen- 
tation of the Llewellyn-Smith formalism using the latest BBA form factors p I] as default. 
The production of baryon resonances, both neutral and charged current, is described us- 
ing an implementation of the Rein-Sehgal model [18] for which 16 baryon resonances are 



4 With the notable exception of coherent pion production, which by its very nature is describing 
scattering off the whole nucleus. 

5 We describe only those of relevance to charged current interactions. 
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included 6 . Coherent pion production is modeled using the Rein-Sehgal model [I i] with 
an updated PCAC formula from a recent revision [15] to the model that takes into ac- 
count lepton mass terms. DIS interactions are calculated using an effective leading order 
model, with modifications at low Q 2 suggested by Bodek and Yang [24]. 

Hadronization is simulated using the AGKY model [25]. It integrates an empirical low- 
invariant mass model with PYTHIA/JETSET at higher invariant mass and is tuned 
primarily to bubble chamber data on hydrogen and deuterium targets. There is a smooth 
transition between the models to ensure continuity of all simulated observables. 

The effect of the nuclear environment is taken into account using an implementation 
of the Fermi Gas model with a modification by Bodek and Ritchie to include nucleon- 
nucleon correlations [26]. Other factors, such as Pauli blocking and the differences between 
nuclear and free nucleon structure functions are also taken into account. Intranuclear 
hadron transport is handled by the INTRANUKE/hA model. It is an effective, data- 
driven, model based on a wide range of hadron-nucleus and hadron-nucleon data. The 
model is validated through comparison to both pion and neutrino scattering data for 
nuclear targets. 

The results of the simulation are shown in Table 1. For the simulation the GENIE default 
values for axial mass of M% EL = 0.99 GeV, and M% ES = 1.12 GeV, were used. The table 
shows the occupancy of primary and final state pion topologies. 





Primary Hadronic System 




Final 
State 


Ott 




7T + 


7T 


2tT 


2tt+ 


27T~ 








> 3vr 


Total 


Ott 


261866 


9187 


38142 





21 


54 





161 





54 


1 


309486 


7T 





32682 


7085 





127 


14 





546 





22 


2 


40478 


7T+ 


549 


890 


139726 





3 


384 





561 





157 


11 


142281 


7T~ 





761 








3 








10 





164 


1 


939 


2n° 














255 








93 





2 


3 


353 


2tt+ 





1 


150 








988 





41 








11 


1191 


2tt- 














1 




















1 


7r°7r+ 


542 


194 


610 





10 


58 





2404 





32 


36 


3886 


7r°7r- 














8 








2 





36 





46 




237 




















48 





594 


33 


912 


> 3tt 





47 


161 





1 


1 





2 








215 


427 


Total 


263194 


43762 


185874 





429 


1499 





3868 





1061 


313 


500000 



Table 1 

Occupancy of primary and final state hadronic systems for a 500,000 event GENIE (release 
2.5.1 was used) sample of on 16 for CC interactions only. A Fermi gas nuclear model and 
a default value of M® L = 0.99 GeV were used. The primary and final state systems were 
separated into different topological groups based on the number of pions. 



6 These are the 16 resonances, of the 18 listed in the original paper, listed as unambiguous in 
the latest PDG baryon tables [23]. 
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3.2 NEUT 



NEUT is able to simulate neutrino interactions from 100 MeV up to a few TeV. The 
Rein-Sehgal model [14], [18] is used to simulate resonance and coherent pion production. 
The GRV94 and GRV98 pdfs [ "] with Bodek-Yang corrections [28] are used to describe 
DIS events. Finally, QE events are simulated using the Llewellyn-Smith [ ] and Smith- 
Mo niz [9] models. There are two recommended default values for the axial mass of 1.11 
GeV or 1.21 GeV; these apply to both QE and RES interactions. For the simulations 
presented here a value of 1.11 GeV was used for both. Nucleon rescattering and me- 
son interactions (especially low momentum pions) in the nucleus are also modeled. To 
describe nuclear effects NEUT uses the cascade model. Each particle is tracked in the 
nucleus until it escapes. For low momentum pions (p < 500 MeV) mean free paths for 
absorption and inelastic scattering are calculated using the L.L Salcedo et al model [29]. 
Higher momentum pion (p > 500 MeV) parameters are taken from experimental results. 
Nucleon re-scattering is simulated considering elastic scattering and single/double pion 
production. 



The results of the simulation using the NEUT generator are presented in Table 2. As 
before the table shows the various pion topologies before and after secondary interactions. 





Primary Hadronic System 




Final 
State 


Ott 




7T + 


7T 


2tt° 


2tt+ 


2tt- 


7T°7r + 


7r°7r- 


7r"*"7r 


> 3tt 


Total 


Ott 


328380 


7148 


33582 





195 


211 





775 





491 


1 


370783 


7T° 


2224 


14498 


7119 





616 


75 





1447 





172 


13 


26164 


7T+ 


2899 


1709 


80376 





60 


702 





1360 





812 


12 


87930 


7T~ 


743 


1632 


742 





66 


6 





126 





899 


9 


4223 


O^O 
Z7T 


11 


99 


92 





711 


5 





270 





13 


15 


1216 


2tt+ 


2 


9 


313 





4 


872 





245 





16 


7 


1468 


2tt- 


2 


6 


2 





6 








2 





6 


1 


25 


ir°ir+ 


14 


114 


508 





115 


141 





3395 





150 


27 


4464 


7T°7r- 


10 


71 


44 





113 


2 





47 





195 


7 


489 


7r~*~7r — 


31 


78 


259 





12 


16 





270 





2195 


21 


2882 


> 3tt 


1 


49 


142 





8 


5 





30 





23 


98 


356 


Total 


334317 


25413 


123179 





1906 


2035 





7967 





4972 


211 


500000 



Table 2 



Occupancy of primary and final state hadronic systems for a 500,000 event NEUT sample of 
on 16 for CC interactions only. A Fermi gas nuclear model and a default value of M(fi EL = 1.11 
GeV were used. The primary and final state systems were separated into different topological 
groups based on the number of pions. 



3.3 FLUKA 



FLUKA simulates the transport and interaction of particles with a focus on the de- 
scription of nuclear models. Particle interactions are described using a Generalized In- 
traNuclear Cascade (GINC) [30]. FLUKA now also describes neutrino interactions. QE 
processes have been included since 1997 and recently the capability to simulate DIS and 
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RES neutrino interactions has been added using the NunDIS and NunRES generators. 
Hadron-nucleon interactions and nuclear effects are taken into account using the main 
reaction mechanism, called PEANUT [30]. For example, pion-nucleon interactions can 
proceed through the non-resonant and the p-wave resonant channel, with the formation 
of a A resonance. The resonance can then either interact with other nucleons or decay. 
These effects can lead to pion absorption or to elastic scattering and charge exchange. 

We present the results of the simulation for FLUKA in Table 3. It is important to mention 
that the recent DIS generator in FLUKA is a beta version and some of the pion production 
numbers in Table 3 are expected to be quite high. Also there is a difference between this 
sample and those produced using the other generators because at present FLUKA does 
not simulate COH pion production and so this process is left out. 





Primary Hadronic System 




Final 
State 


Ott 




7T + 


7T 


2n° 


2tt+ 


2tt- 


7T°7T + 


7r°7r- 


7T~*~7T 


> 3tt 


Total 


Ott 


267108 


13702 


41299 





47 


4 





129 





25 





322314 


7r° 


964 


38871 


8453 





307 


4 





465 





9 





49073 


7T+ 


2005 


2903 


116159 





16 


29 





421 





91 


2 


121626 


7T _ 


121 


1989 


438 





18 








26 





89 





2681 


277° 





115 


122 





548 


1 





68 











854 


2tt+ 


2 


11 


230 





2 


37 





61 








2 


345 


2tt- 





1 


1 


























2 


7r°7T + 


2 


177 


391 





53 


6 





1383 





21 





2033 


7T°7T- 





47 


15 





38 








9 





14 





124 


7T+7T~ 





118 


320 





1 








78 





392 


1 


910 


> 3tt 





7 


13 





1 








5 








12 


38 


Total 


270202 


57941 


167441 





1031 


81 





2645 





641 


18 


500000 



Table 3 

Occupancy of primary and final state hadronic systems for a 500,000 event FLUKA sample of 
Ua on 16 for CC interactions only. The primary and final state systems were separated into 
different topological groups based on the number of pions. 



3.4 NuWro 

NuWro uses the Llewellyn Smith model [ I 1 ] with the latest BBA form factors [ I] for QES 
events. DIS processes are described using the GRV94 pdf with Bodek-Yang corrections. 
The Rein-Sehgal models [14] and [15] are used to describe COH pion production. We now 
give a more detailed description of how these processes are implemented in NuWro [31] : 

• The RES region is defined by a cut at W < 1.4 GeV, where W is the invariant hadronic 
mass and has the standard definition 7 . Only A resonance is considered with form 
factors from a fit to ANL and BNL single pion production data [32] . The non-resonant 
background is modelled as a fraction of the DIS contribution. 



7 W 2 = (P + q) 2 , where P is initial nucleon 4-momentum and — q 2 = Q 2 is the squared 4- 
momentum transfer. 
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• In the region of 1.4 GeV < W < 1.6 GeV, the A contribution is linearly turned off as 
the DIS contribution is turned on. 

• The DIS models are applied when W > 1.6 GeV. Together with the Bodek-Yang model 
NuWro's own hadronization model is used. 

• The COH pion production contribution, as predicted by the Rein-Sehgal model, is 
multiplied by a factor of 0.66. 

NuWro offers a choice of two basic nuclear models: a relativistic fermi gas (RFG) model 
or an effective spectral function model. For the tables shown here the RFG model was 
used. We also used the default value for the axial mass of Ma = 1.1 GeV. In the appendix 
Tables 8 and 9 show results for samples generated with a non-default axial mass and for 
both the RFG and effective spectral function nuclear models. 

3.4-1 NuWro's Own Intranuclear Cascade Model 

For this report two different intranuclear cascades were used, one included with NuWro 
and another, the Bertini cascade, that comes from Geant4. First we show the results 
using NuWro's intranuclear cascade in Table 4. 





Primary Hadronic System 




Final 
State 


Ott 




7T + 


7T 


2tT 


27T+ 


2vr- 


7r°7T + 


7r°7T- 


TPTT 


> 3?r 


Total 


Ott 


281781 


10437 


49381 


1 


49 


78 





523 





223 


2 


342575 


7T° 


2137 


17826 


12404 





178 


29 





1095 





85 


24 


33778 


7T+ 


5241 


3168 


97356 





32 


255 





1187 





407 


4 


107650 


7T — 


1165 


2849 


1667 


1 


22 


4 





140 





413 


4 


6265 


2tt° 


12 


332 


123 





179 


4 





203 





13 


38 


904 


2tt+ 


11 


39 


1948 





3 


248 





245 





10 


4 


2508 


2tt- 





20 


2 





1 








1 





11 


2 


37 


7T°7T+ 


26 


382 


1090 





38 


65 





1963 





100 


18 


3682 


n°n- 


7 


237 


94 





37 








55 





76 


7 


513 




22 


67 


783 





6 


8 





220 





804 


8 


1918 


> 3tt 


1 


27 


47 





10 


11 





69 





29 


76 


270 


Total 


290403 


35384 


164895 


2 


555 


702 





5701 





2171 


187 


500000 



Table 4 

Occupancy of primary and final state hadronic systems for a 500,000 event NuWro sample of 
on 16 for CC interactions only. Sample generated with NuWro vll2. M^ ELcc = 1.10 GeV, 
nuclear model: Fermi gas. The primary and final state systems were separated into different 
topological groups based on the number of pions. 

3.4-2 Combination of NuWro and Geant4 generators 

Now we show the results for the simulation using the Geant4 [33] intranuclear cascade. 
The primary neutrino interactions on oxygen are generated in NuWro and secondary 
particles are propagated inside nucleus using the Bertini cascade [34]. 

There are two intranuclear cascades in Geant4: Bertini and Binary. However we only 
consider the Bertini cascade. It works very well for particles with energy below 3 GeV. It 
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takes into account a variety of interactions that nucleons and pions can undergo, using 
very recent cross-sections. Compared with Binary, the Bertini cascade is better tested and 
developed by both members and users of Geant4. A detailed description can be found 
in [34]. 

We would like to state that, in its original form, the Bertini cascade was not set up 
to propagate secondary particles from neutrino interactions. A particle from a neutrino 
interaction can originate at any point within the volume of the nucleus. This meant 
changing the original Bertini cascade, in which a particle always hits a nucleus from the 
outside (so that the cascade always starts on the surface of the nucleus) so that the point 
where the particle enters is selected uniformly over the volume of the nucleus. The results 
of the simulation using this approach are shown in Table 5. 





Primary Hadronic System 




Final 
State 


Ott 




7T + 


7T 


27T 


2tt+ 


2vr- 


7r°7T + 


7r°7T- 


7P7T 


> 3?r 


Total 


Ott 


282443 


8888 


36738 





67 


47 





616 





173 


18 


328990 


7T° 


3148 


20328 


8963 





147 


25 





801 





65 


25 


33502 


7T+ 


3866 


2770 


115609 





31 


218 





1162 





323 


6 


123985 


7T — 


908 


2845 


1190 


2 


34 


1 





134 





430 


10 


5554 


Z7T 


8 


123 


152 





208 


1 





91 





4 


24 


611 


2tt+ 


7 


32 


676 





1 


363 





199 





1 


1 


1280 


2tt- 





8 


3 





4 








3 





9 





27 


7T°7r+ 


12 


161 


904 





30 


41 





2448 





53 


10 


3659 


TT°TT~ 


2 


89 


70 





26 








30 





49 


9 


275 




9 


135 


578 





5 


4 





197 





1056 


22 


2006 


> 3tt 





5 


12 





2 


2 





20 





8 


62 


111 


Total 


290403 


35384 


164895 


2 


555 


702 





5701 





2171 


187 


500000 



Table 5 

Occupancy of primary and final state hadronic systems for a 500,000 event NuWro (plus 
GEANT4) sample of Vn on 16 for CC interactions only. Sample generated with NuWro vll2. 
M® ELcc = 1.10 GeV, nuclear model: Fermi gas. Also, hadrons produced inside the nucleus 
were propagated using the GEANT4 cascade model. The primary and final state systems were 
separated into different topological groups based on the number of pions. 
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4 Discussion and Conclusions 

Looking at the tables for each generator we can see some large differences in both the 
primary and final state topologies (number of pions at the initial and final states). In 
many cases these differences are above statistical fluctuations 8 and reflect the difficulty 
in modeling this energy region. At ~ 1 GeV we are in a transition region where QE, RES 
processes dominate but where there is also a significant DIS component being switched on 
as we increase in energy. Although the models used to describe these processes separately 
are often common to a generator, there are still many differences in the way in which a 
particular generator handles the merging of the relative contributions in this transition 
region. This combined with differences in the assumed nominal values for many of the 
input parameters can lead to different predictions from the same set of models. 

In general QE processes give rise to topologies with no pions in the initial and final states 
whereas DIS and RES are more likely to result in events with pions in the primary and 
final state. All the generators have a larger number of On topologies in the final state 
than were in the primary state. This indicates that pions are more likely to be absorbed 
than created. One noticeable feature is that NEUT has a larger number of On topologies 
in both the initial and final states than the other generators. This is an indicator that 
either the cross-section for QE processes in NEUT is higher than in the other generators 
or that the contribution from DIS and RES events is lower. Mainly we expect differences 
to depend on the values of form-factor parameters and the nuclear models used in each 
generator 9 . 

Resonant events are tricky to handle; there are large uncertainties on the underlying 
cross-sections as well as differences in the way each generator models the propagation of 
pions in the nuclear environment. Which model is used is important and we can see the 
effect of these differences in the tables presented in this report. For instance, if we look 
at the On primary state column we see that the number of these events which result in 
a single n + in the final state topology varies from 549 (Table 1) all the way up to 5241 
(Table 4). Still looking at the first column we see that the number of 7r°'s and 7r~'s in the 
final state can vary from zero (Table 1) to 1-3 thousands (Tables 2, 3, 4, 5). All of these 
events had no pions in the primary state and so these numbers directly reflect differences 
in the hadron transport models used in the generators. 

As mentioned before the tables presented so far contain important information about 
FSIs. To elucidate the effects of FSIs we have compiled a summary table (Table 6) showing 
directly the topology changing effect due to intranuclear hadron transport. For each of 
the generators, and for a selection of primary and final state topologies, we show, out 
of all events with a given primary state topology, the fraction of those which have both 
the primary and final state topology. Looking at Table 6 the first two rows tell us the 
fraction of events, with a single pion in the initial state, that will still have a single pion 
in the final state. We can see that the most probable scenario is that a pion created at the 
primary vertex will not re- interact (we can look at this as a transparency of the nucleus) . 



8 In some of the cells with enough high statistics these differences can be above 10%, well above 
statistical fluctuations. 

9 Or, if available, the choice of simulation within a particular generator. 
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The next two rows show the percent of pions which are absorbed and the remaining 
rows show the effect due to charge exchange processes (we assume that a pion, to first 
order, re-scatters only once). It seems that NuWro has a lower transparency compared to 
the other generators, whilst GENIE's is higher than that of the others. This may be in 
response to the absorption and charge exchange processes, for which GENIE could have 
too little and NuWro too much. Despite these differences the agreement is still good, 
considering the complexity of FSIs, and one can convince oneself that the analyzed MCs 
give quite similar results. These modes are important to neutrino oscillation experiments 
since single pion events form the main background channel. 





GENIE 


NEUT 


FLUKA 


NuWro 


NuWro + G4 




75% 


57% 


67% 


50% 


57% 


TT + — > 7T + 


75% 


65% 


69% 


59% 


70% 


7T° -> Ovr's 


20% 


28% 


24% 


29% 


25% 


7r + — * Oir's 


20% 


27% 


25% 


30% 


22% 




2% 


7% 


5% 


9% 


8% 




2% 


6% 


3% 


8% 


8% 




4% 


6% 


5% 


8% 


5% 



Table 6 

Rate of events with single pion or no pion in final state if there was single pion in initial state. 

Taking into accout recent MiniBooNE results [35] of CCl7r + to CCQE cross section ratios 
we can approximately compare the corresponding numbers from our simulations with each 
other and with the MiniBooNE measurements. We present this comparison in Table 11 
in the appendix. 

In this report we have presented event numbers, with an emphasis on pion production, 
for a large number of the current neutrino generators. We have started to discuss the 
differences between the different generators and models. These differences are now a 
matter for further investigation. A more extensive set of tables has been published on 
the Ladek MC Project website and more tables, at different energies, are planned to 
be produced in the near future. We hope that the tables will be a valuable resource 
to the reader. Improvements in the theoretical models used by generators, and better 
measurements made by future neutrino experiments, will, by improving our understanding 
of neutrino cross-sections, help shed light on the elusive nature of the neutrino. 
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5 Appendix 

In this appendix we present some tables showing the effect of changing the values for the 
axial mass parameter or the use of a different nuclear model. We also present some results 
from the Nuance generator which, due to problems with generating samples comparable 
to those produced using the other generators, were left out of the main report. Finally we 
show some comparisons to a MiniBooNE measurement of the CC single pion-production 
to quasi-elastic cross-section ratio. 

GENIE 
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Final 
State 
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7T + 
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2tt- 




7r°7T- 
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6 





1 
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Total 


287741 


39869 


165663 





379 


1456 





3591 
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312 


500000 



Table 7 

Occupancy of primary and final state hadronic systems for a 500,000 event GENIE sample 
of i/„ on 16 O for CC interactions only. A Fermi gas nuclear model and a non default value 
of MaQEL = 1.18 GeV were used. The primary and final state systems were separated into 
different topological groups based on the number of pions. Coherent pion production events 
were counted as having a single pion in the final state only. 

NuWro 

There are two additional tables for NuWro event samples. In Table 8 a Fermi gas nuclear 
model was used but with an axial mass equal to 1.0 GeV. In Table 9 an effective spectral 
function was used instead of the Fermi gas model. 

Nuance 

Nuance is an advanced and freely available neutrino generator written by Dave Casper 
of the University of California (the version used in this study is 3.006). The program is 
most suitable for generating events on oxygen/water, as its FSI model was thoroughly 
tested for this target. Resonant and coherent/ diffract ive interactions are simulated using 
the Rein-Sehgal model. More details can be found in [7]. Unfortunately in the case of 
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Table 8 

Occupancy of primary and final state hadronic systems for a 500,000 event NuWro sample of 
v„ on 16 for CC interactions only Sample generated with NuWro vll2. M® ELcc = 1.0 GeV, 
nuclear model: Fermi gas. The primary and final state systems were separated into different 
topological groups based on the number of pions. 
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Table 9 



Occupancy of primary and final state hadronic systems for a 500,000 event NuWro sample of 
f M on 16 for CC interactions only. Sample generated with NuWro vll2. M® ELcc = 1.0 GeV, 
nuclear model: effective spectral function. The primary and final state systems were separated 
into different topological groups based on the number of pions. 

Nuance it was not possible to create a single sample consisting of events with final state 
particles before and after the FSIs have been applied. Therefore we were unable to produce 
tables like those presented for the other generators. Instead, in Table 10 we present pion 
statistics for two separate samples independently, where one had FSIs turned off and in 
the other they were left on. 
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Table 10 

Two samples, for muon neutrino interactions on oxygen, of 500,000 events each were used, one 
with FSIs turned on and one with FSIs turned off. The muon neutrinos had an initial energy 
of 1 GeV and only CC interactions were taken into account. A Fermi gas nuclear model was 
used. The default axial masses used were Ma = 1.1 GeV for resonant single pion production, 
Ma = 1.3 GeV for resonant multi-pion production and Ma = 1-03 GeV for coherent/diffractive 
pion production (the same value was used in QE event generation; the vector mass for QE was 
My = 0.84 GeV). The samples consist of QE(59%), RES(39%) and COH(2%) events. 

The single charged pion production to quasi-elastic cross section ratios 

In Table 11 the comparison between the ratio of l7r + /07r is shown for the generators we 
considered. From the MiniBooNE data [35] re-scaled for an isoscalar target and corrected 
for FSI's (meaning events at the initial vertex and before any hadronic re-interactions) 
the following values of the ratio are 10 : 0.52 ± 0.06 for a neutrino energy of 0.95 ± 0.05 
GeV and 0.63 ± 0.07 for a neutrino energy of 1.05 ± 0.05 GeV. These numbers agree well 
with the corresponding numbers, see the PHS total column in Table 11, obtained in our 
simulations, even though the Genie ratio is slightly higher and the Neut ratio is smaller. 



Our simulations, Charged Current interactions on 16 


E v = 1 GeV 


l7r+/07r (PHS total) 


17T+/07T (Final State total) 


GENIE 


0.706 


0.460 


NEUT 


0.368 


0.237 


FLUKA 


0.620 


0.377 


NuWro 


0.568 


0.314 


NuWro+G4 


0.568 


0.377 


Nuance 


0.579 


0.359 



Table 11 

The ratios l7r + /07r obtained in our simulations. 



The MiniBooNE ratios of CClvr+/CCQE are taken from Fig. 2 [35]. 
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